Topological phases and delocalization of quantum walks in random environments 
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We investigate one-dimensional (ID) discrete time quantum walks (QWs) with spatially or tem- 
porally random defects as a consequence of interactions with random environments. We focus on the 
QWs with chiral symmetry in a topological phase, and reveal that chiral symmetry together with 
bipartite nature of the QWs brings about intriguing behaviors such as coexistence of topologically 
protected edge states at zero energy and Anderson transitions in the ID chiral class at non-zero 
energy in their dynamics. Contrary to the previous studies, therefore, the spatially disordered QWs 
can avoid complete localization due to the Anderson transition. It is further confirmed that the 
edge states are robust for spatial disorder but not for temporal disorder. 



PACS numbers: 03.67.-a, 05.40.Fb, 73.20.Fz, 03.65.Yz 



I. INTRODUCTION 

A quantum walk (QW)0 describes quantum mechan- 
ical time evolution of particles, which is identified as a 
random walk when the system is brought to the classical 
limit. The QW may provide a unique avenue to realize 
quantum computation since the QW can be applied for 
efficient algorithms of quantum computationJa^EI Among 
several kinds of QWs, the two-state discrete-time QW in 
one dimension (ID) has been intensively investigated due 
to its simple formalism. More remarkably, it has been 
experkpentally realizad-iin various systenie, such as cold 
atomga, trapped ionsBfl, and photons.ErB The discrete- 
time QW is described by two basic operators. The coin 
operator C defined as 



C = 



cos 9 
sin 6* 



— smt 
cos 9 



(1) 



acts on two internal states, right and left walkers \R) and 
\L), respectively. The shift operator S is defined as 



N/2-1 



{\n+l){n\(S)\R){R\ + \n-l){n\(^\L){L\), (2) 



n=-N/2 

SO that each walker moves from a position n to its corre- 
sponding neighbor n ± 1 . Here, N represents the system 
length. The time evolution operator U is built up from 
these two operators and the corresponding Hamiltonian 
H is defined through 



U = S(j:jn){n\<E>C) =e 



iHSt 



(3) 



where St is a unit of time and we hereafter set 6t = 
1. Thus, a state at time t, \ip{t)), is given by \i^{t)) = 
U*\ifj{Q)) starting from an initial state \ip{0)). 

It is known that the Hamiltonian H described by Eqs. 
(0) " (H) iS[eaiuvalent to the Dirac equation in the contin- 
uum limitli3'U Recently, symmetries of the Haniiltonian, 
which are relevant to classify topological phaseaiSllj, are 
examined in Rcf. and it is found that the QW de- 
scribed by Eqs. (0)-(P), which has been realized in many 



experimentJi^i, possesses chiral symmetry. Ref. clar- 
ified that the Dirac equation derived from the ID QW 
with chiral symmetry gives a finite Berry phase and 
thus generates edge states near boundaries, since ID 
chiral classes caiiJpe characterized by an integer topo- 
logical numberEjtJ. The edge states in the topologi- 
cal phase wrt«Jd-|be also useful for topological quantum 
computationll§tO, if the system is described by Majorana 
fermions, which is currently a subject of intensive studies 
in condensed matter physics. 

Understanding effects of the interaction with environ- 
ments giving rise to decohererice is a key issue to real- 
ize the qupiiiim computation£3_jErLthis end, there are 
theoreticaltfrO and experimentalmLl studies of discrete- 
time QWs taking account of spatial and/or temporal dis- 
orders. It is found that Anderson localizations occur for 
QWs with spatial disorder while QWs with temporal dis- 
order approach to classical random walks. 

In this work, we investigate two-state discrete time ID 
QWs described by Eqs.(^-(||), which belong to the chiral 
orthogonal class, interacting with spatially or temporally 
random defects. Contrary to the previous studies, we re- 
veal that the QWs with spatial disorder exhibit delocal- 
ized behaviors. This remarkable conclusion is drawn by 
extensive numerical calculations, and is further ensured 
by symmetry arguments for the ID chiral class. We also 
find that the edge states in the topological phase are ro- 
bust for spatial disorder but not for temporal disorder. 

This paper is organized as follows. In the Sec. II, we 
give general remarks on QWs without any disorder and 
demonstrate that the edge states are induced by intro- 
ducing a refiecting coin operator. Our main results for 
the spatially or temporally disordered QWs are presented 
in Sec. HI. Section IV is devoted to discussions and sum- 
mary. 



II. QUANTUM WALKS WITH REFLECTING 
COIN OPERATORS 



Let us first explain eigenstates of the QW. Eigenenergy 
uj\ of the QW, which has periodicity of 27r, is defined 
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FIG. 1: (Color online) The QW with a reflecting coin oper- 
ator at the origin. 



through the eigenvalue of the time-evolution operator by 
U\ip\) = e*"^^ |'0a), where \tp\) is the corresponding eigen- 
state. Hereafter, the subscript A is written only when 
we need to emphasize. The dispersion relation result- 
ing from Eqs. (|l|)-(||) is given_hjL cos(w) = cos(fc) cos(0), 
where fc is a wave number.Eilo Since the energy gap 
should exist to support edge states in topological phases, 
edge states of the QW are able to appear when 9 ^ 0,tt. 

To generate the edge states, a boundary should be pre- 
pared for the QW. While a split-step method is proposed 
in Ref. ^ we employ a simpler method to make the 
boundary in terms of a reflecting coin operator defined 
as 



Tl 
±1 



(4) 



so that a right walker is changed to a left walker and 
vice versa. Thereby, realizes a hard wall boundary. 
The sign of the reflecting coin oper|a±aL determines the 
presence or absence of edge states £jo In case of the 
Hadamard walks {6 — 7r/4), the edge states appear if Cj^ 
is introduced at the boundary. 

We now consider the Hadamard walks with Cj^ at the 
origin as shown in Fig. ^ The initial state is prepared as 

\m)^\o)(s{\R)+m)/V2 (5) 

throughout the paper, so that the right and left walkers 
do not mix with each other due to real-number elements 
of the coin operator and then the probability distribution 
becomes symmetric in the absence of spatial disorder. 
Figure || (a) shows the probability distribution Pn{t) at 
t = 80, which is defined as 



a=R,L 



(6) 



We find a sharp peak at the origin due to edge states in 
the topological phase and two smaller peaks near n ^ 
±50 commonly observed in QWs.l3a. Note that these 
observations are consistent with the recently derived limit 
measurctZl at i — > od for the ID QW on a half line. The 
time dependence of the survival probability 



Po{t) = Pn{t)\n=0 

and the position variance of QWs 

v{t) ^ {m\n'\m) - [{m\n\m)? 



(7) 



(8) 





t 



FIG. 2: (Color online) Probability distributions at t = 80 of 
the QWs with 59s = (a) and 7r/4 (b) with for a single 
sample, t dependence of the survival probability Po{t) (c) and 
the variance v{t) (d). In (c) and (d), the solid and dashed 
curves represent the QWs with 58s ~ and 7r/4, respectively. 
The two thin lines in (d) indicate power law behaviors. The 
number of samples is 10* in case of disordered QWs. 



are shown by solid curves in Fig. ^ (c) and (d), respec- 
tively. We confirm that the edge states are robust for 
time evolution since Pq {t) keeps its value unchanged and 
v{t) is,jproportional to t^. which agrees with ordinary 
QWs.ffl 



III. DISORDERED QWS 

Next, we consider disordered QWs. We assume that 
disorder is introduced via fluctuations of 9 in Eq. (|^), 
which do not break chiral symmetry. To this end, we 
redefine 9 so that 9 is randomly distributed over positions 
n or time t as 

ee [9~69sit)/'2:9 + 69,(^t)m, 

where 9 represents the mean value of distributed 9 and 
59s{t) indicates the strength of spatial (temporal) disor- 
der. We hereafter restrict our attention to QWs with 
9 ~ tt/A and with cither spatial or temporal disorder. 



A. spatially disordered QWs 

Let us begin with the spatial disorder. The probabil- 
ity distribution Pn(t) for spatially disordered QWs with 
after 80 time steps is shown in Fig. ^ (b). The peak 
at the origin due to the edge states seems to be robust 
for the spatial disorder. The survival probability Po{t) 
indeed ensures this point: the profile of Po{t) is almost 
unchanged from that for the clean QW. On the other 
hand, the outer two smaller peaks existing in Fig. ^ (a) 
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disappear in Fig. ^ (b), implying that the QWs suffer 
from the spatial disorder except their edge states. In 
this case, the Anderson localizatioji wojild be important 
as studied in the previous works. BE2lLil If the Anderson 
localization completely dominates dynamics of QWs, the 
variance v(t) should become a constant after many time 
steps. However, v{t) for the spatially disordered QWs 
[Fig. H (d)] exhibits a power law behavior with a smaller 
exponent, suggesting that the QWs show anomalous dif- 
fusion even in the presence of spatial disorder. 

To address the above point, we look into the density 
of states (DOS), p{uj), which provides us with a clear- 
cut interpretation of the seemingly peculiar dynamics of 
the QWs. Note that a state Ifpit)) is considered as the 
superposition of all the eigenstates because we employ a 
delta-function like initial state \ip{0)). The DOS for the 
clean QW with translational symmetry is given as 



27r\/cos2 6 — cos^ ( 



for cos^ uj < cos^ 9. (9) 



When cos^ uj = cos^ 9 at 9 ^ 0, tt, pd^j) diverges due to 
the van Hove singularity. For the QW with , we find 
that the DOS consists of Pc(w) and doubly degenerate 
states at oj = 0, TT as shown in Fig. ^ (a). According to 
the nature of for uj\ ~ 0, tt, these states are identified 
as the edge states in the topological phase. In the clean 
QW, therefore, the edge states appear not only at cj = 
but also at w = TT. With increasing S9s [Fig. ||(b)], the 
divergences in the DOS at w = ±7r/4 and ±37r/4 are 
rounded since the van Hove singularity is not well defined 
for disordered systems. Furthermore, we notice that the 
energy gaps around w = 0, tt are reduced while the delta 
functions due to edge states robustly remain. 

Remarkably, there appear additional divergences in the 
DOS at w = ±7r/2 for spatially disordered QWs [Fig. | 
(b)]. We give attention to the divergences in p{uj) at 
Lo — ±7r/2 by showing the enlarged picture of DOS near 
Lo = 7r/2 in Fig. ^ (d). We find that the divergence be- 
comes stronger with increasing 69^ and the DOS for the 
smallest disorder (69 s ~ tt/S) shows the oscillating be- 
havior. This oscillation reminds us p£ the DOS derived 
from chiral random matrix theories. Indeed, a system 
belonging to the ID chiral class is known to show an 
Anderson transition at zero energyr-^hich is accompa- 
nied by the divergence in the D0S.E3 Although the en- 
ergy at which the DOS diverges is different from what 
we observed for the QW, we claim that the Anderson 
transition in the ID chiral class occurs in the spatially 
disordered QW at lo = ±7r/2. This is supported by the 
following symmetry arguments. Chiral symmetry gives 
a constraint that eigenenergies appear in pairs of oppo- 
site sign and this makes the zero energy states in chiral 
classes singular. In the QWs, chiral symmetry leads to 
the eigenvalues of U appearing in pairs like e^^'^K In 
addition, since S in Eq. gives only nearest neighbor 
hopping, U in Eq. has the bipartite structuije and the 
eigenvalues of U also appear in pairs as ±e*"* .EJ Thereby, 
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FIG. 3: (Color online) DOS p{uj) for the QWs with for 
SBs = (a), SOs = 7r/8,7r/4,7r/2, and tt (b), and 2tv (c). 
N = 500 and the number of samples is 10^ for each SOs{^ 0). 
(d) The enlargement of DOS around u = tt/2 for QWs with 
various S6b. (e) The scaling collapses into Eq. (hol) for the 
DOS near uj = tt/2 for QWs of the larger system A'^ — 10^ 
with IC* samples. The labels in (d) represent values of 50s. 
In calculations of DOS, the periodic boundary conditions are 
imposed on both edges of the ID QW. 



there are four eigenenergies related to each other. 

Redefining ujx as wa = ^\ + t^/'2, these eigenenergies are 
rewritten as 



{cj'^ + tt/2, -uj'^ - tt/2, uj'^ + 3^/2, 



tt/2). 



One finds that the first and fourth (second and third) 
eigenenejigies appear in pairs with respect to w = tt/2 
(— 7r/2)Eil, making lo ~ ±tt/2 special energies as in zero 
energy of ordinary chiral classes. We therefore come to 
the conclusion that the Anderson transition of the ID 
chiral class should occur at a; = zLtt/2 even if there are 
energy gaps and the edge states exist at w = 0, tt. 

According to universality, the same critical behaviors 
near zero energy studied with a tight-binding Hamilto- 
nian in Ref. g9| should be observed in the discrete time 
QW near lo = zLtt/2. Thereby, we write down the mod- 
ified formulae for divergences of the DOS p and the lo- 
calization length ^, defined as the inverse of Lyapunov 
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FIG. 4: (Color online) (a) Dependence on Suj of the localiza- 
tion length ^ of the QWs with 59s — 7r/8, 7r/4, n/2, n, and 27r 
from the top to the bottom, (b) The scaling collapse into Eq. 
(^) for ^ in (a). The parameters and r are obtained by 
the fitting. The dotted thin line represents Eq. (0). 

exponents, at w = ±7r/2 as 

p{io) = po{6u;T\\n\6u;T)\y\ (10) 
e(w) - Co|ln((5c.r)|, (11) 

respectively, where du = |a;^7r/2| represents the distance 
from the critical points uj = ±7r/2. r denotes the mean 
free time and po and are constants. As shown in Fig. 
0(e), we verify that the DOS near w ~ 7r/2 obeys Eq. 
(|lO|) fairly well. 

We also calculate the localization length ^ of the spa- 
tially disordered (JW near lj = 7r/2 by using the transfer 
matrix mcthodE3'E3 By rearranging the wave function 
amplitudes around the position n of eigenvalue equation 
U\ip) = e*'^|'0) , we obtain the iterative relation, 

rp _ { e"*"/cos6'„ -tan6'„ \ , , 
" ~ -tan6'„ e'^^/cosOn)' ^ "^^ 

where '>pn,iT=R.L represents the wave function amplitude 
of the right or left walker at the position n and 0„ denotes 
6 at the position n. The huge system length N — 10^ 
allows us to evaluate accurate numerical values of ^. 

Figure ^ (a) shows the Suj dependence of ^ for various 
66s- We confirm that ^ for any 69 s shows the diverging 
behavior in the vicinity of (5a; = while ^ for 7^ is 
finite, resulting in the Anderson localization. Although 
the numerical restriction of 6lu > 10~^^ prevents us to 



reach £ — > c» with accuracy, we confirm that ^ follows 
Eq. ( |ll|) as shown in Fig. ^ (b). This scaling collapse gives 
an evidence that diverges at 6uj = for any 66s- In this 
case, the transport property is anomalous by involving al- 
gebraic decays instead of the exponential ones.EJ Hence, 
the walkers can move around through these delocalized 
states at cj = ±7r/2 and then v{t) indicates anomalous 
diffusion as shown in Fig. ^ (d). Interestingly, three kinds 
of states (edge, localized, and delocalized states) are si- 
multaneously appear in dynamics of spatially disordered 
QWs with Cj^- This coexistence occurs not only for the 
specific initial state, Eq. (|^), but also for general delta 
function like initial states, so that dynamics of QWs is 
described by the whole eigen states. We mention that 
since the coin operator of the, previous works for the 
ID spatially disordered QWaJ'E2rEll does not retain chiral 
symmetry, the Anderson transition is not found and only 
the Anderson localized states arc confirmed. 

It is known for zero energy states of chiral classeJsl and 
various topological insulators, such as the integer quan- 
tum Hall transition, that topologically protected edge 
states stably exist unless the bulk band gap is closed. 
When a symmetry-preserving perturbation collapses the 
bulk gap, the delocalized states which can be coupled 
with the edge states so that the topological number is 
changed should appear even in ID disordered systems. 
In other words, the existence of the delocalized states 
as the Anderson delocalization is essential to topological 
phases in disordered systems. It is naturally expected, 
then, that the ID QW with chiral symmetry also follows 
a similar fashion under the spatial disorder of 6. 

Here, we clarify how the strong spatial disorder 66s 
induces the delocalized states at w — 0,7r as a conse- 
quence of the gap closing. Figure ^ (b) indicates that 
the bulk gaps around uj — O^n become smaller with in- 
creasing 66 s, though there still remain the edge states 
at 66s = TT. Unfortunately, it is difficult to find the gap 
closing from the DOS since the finite system calculation 
keeps the mean level spacing finite. Alternatively, we 
calculate the localization length ^ at = for the QW 
without Cj^ for various 66s by using Eqs. (|l2|) and (|l|). 
Note that when the bulk gap around a; = is finite, the 
localization length is interpreted as of an exponentially 
decaying evanescent mode. Figure ^ clearly shows that 
the localization length ^ at w = increases with increas- 
ing 66s and finally diverges at 66s = 27r, where 6 is fully 
distributed in 2tt radians. This suggests that the bulk 
gap around w = is closed. Figure y (c) shows the DOS 
p(uj) at 66s = 2tt- We confirm that, instead of the delta 
frmctions of the edge states, the divergences in the DOS, 
which are the same with those at w = ±7r/2, appear even 
at cj = 0, TT. Therefore, the Anderson transitions occur 
at a; = 0,7r as well as ±7r/2 at 66s = 27r and the coexis- 
tence of the edge states and the delocalized states docs 
not occur. 
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(a) 10° 



FIG. 5: (Color online) SOs dependence of the inverse of the 
localization length at a; = 0. The length of the system is 
A*' — 10*. Inset: The enlargement of near S9s = 2-k. 



B. temporally disordered QWs 

Finally, we consider QWs with temporal disorder 
where 9 in Eq. (|l]) is shuffled in each time step. It is 
known that, since the temporal disorder gives rise to de- 
coherence, the QWs approach to the classical random 
walks where the probability distribution of QWs becomes 
aian and the position variance u(t) is proportional to 
1 We study how the edge states of the QWs with Cj^ 
are affected by the temporal disorder. The probability 
distribution P„(t) after 80 time steps in Fig. || (a) demon- 
strates that the peaks due to the edge states are substan- 
tially reduced with increasing S9t. For longer time steps 
t ~ lO'*, the sample averaged P„(i) with temporal dis- 
order [Fig. 1^ (b)] become Gaussian, indicating that the 
QWs approach toward the classical random walks. The 
tiny peaks at the origin, which would be the remnants of 
the edge states, arc found for the QWs with C]^ . In Fig. 
(a), the survival probabilities Po{t) of QWs with are 
compared with those without C^. We find that Po{t) of 
the former QWs behaves differently from the latter dur- 
ing rather short time steps, while at longer time steps 
they converge to the same curve, which means that the 
edge states gradually disappear. The position variance 
v(t) in Fig. (b) increasing linearly in time for longer 
time steps clarifies that the QWs with edge states are 
also transformed into the classical random walks. 



IV. DISCUSSION AND SUMMARY 

In this work, we have focused only on discrete-time 
QWs. Here some comments are in order on the rele- 
vance of our results to continuous-time QWs with spa- 
tial disorder. While the ID continuous-time QiiLwith 
on-site random potentials is studied previouslyEjo, the 
corresponding Hamiltonian does not retain chiral symme- 
try. We note that the Hamiltonian of the ID continuous- 
time QW possessing chiral symmetry is the tight-binding 
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FIG. 6: (Color online) Probability distributions at t = 80 for 
a single sample (a) and sample averaged probability distribu- 
tions at t = 10* (b) of the temporally disordered QWs. The 
solid and dashed curves represent QWs with 59t — tv/S and 
7r/4, respectively, and the thick and thin curves distinguish 
the QWs with and without C'^ , respectively. The number of 
samples is 10*. 



model with only nearest neighbor random hopping terms. 
Indeed, the latter model is consistent with the one stud- 
ied in Rcf. ^ and the Anderson transition is possible 
only at zero energy in contrast to the discrete-time QWs 
studied in the present work. Therefore, the coexistence 
of edge states and the Anderson transition is peculiar to 
discrete-time QWs. 

We also mention that, although we have been con- 
cerned with the QWs with here, our conclusion can 
directly be applied for the QWs without , except for 
the argument on edge states at a; = 0, tt. 

In summary, wc have investigated two-state discrete 
time QWs belonging to the ID chiral orthogonal class in 
the presence of spatial or temporal disorder. We have elu- 
cidated that Anderson transitions in the ID chiral classes 
occur at cl) = ±7r/2 in the QWs with any strength of spa- 
tial disorder and thereby the QWs can avoid the complete 
localization. We note that, while the currently available 
experiments on discrete-time QWs can accomplish only 
a few ten time steps and it may difficult to completely 
eliminate the temporal disorder giving rise to decoher- 
ence in the spatially disordered QWs, the delocalization 
behavior found in the present work will be observed when 
the technology of QWs will be further developed in the 
application of quantum computer. 

Furthermore, we have shown that the coexistence of 
edge, localized, and dclocalized states is realized in the 
time evolution of the spatially disordered QW with C]^ . 
This characteristic nature of the QWs is supported by 
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FIG. 7: (Color online) t dependence of the survival probabil- 
ity Po{t) (a) and the variances v{t) (b). The solid and dashed 
curves represent QWs with S9t = n/8 and 7r/4, respectively, 
and the thick and thin curves distinguish the QWs with and 
without , respectively. The number of samples is 10". The 
dotted line in (b) indicates v{t) oc t. 



chiral symmetry of H and bipartite structures of U. We 
predict that delocalization of QWs with spatial dipordcr 
should be observed in a wide variety of ID QWaia be- 
longing to the other chiral classes and class D and Dili 
which are described by Majorana fermions, since these 
uniKecsality classes also show the Anderson transition in 
IDO'E^ The QWs for which tuning system parameters 
and observing walkers' probabilities are possible in ex- 
periments would provide an intriguing arena to study 
topological phases for systems with defects and decoher- 
ence. 

Note: After submittin|g-|the manuscript, we found a 
preprint arXiv:1105.5334e3 in which the edge states of 
QWs are observed in experiments. 
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